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We present and analyze long-term optical photometric measurements of the three active stars V2253 Oph, IT Com and 
IS Vir. All three systems are single-lined spectroscopic binaries with an early K giant as primary component but in different 
stages of orbital-rotational synchronization. Our photometry is supplemented by 2MASS and WISE near-IR and mid-IR 
magnitudes and then used to obtain more accurate effective temperatures and extinctions. For V2253 Oph and IT Com, 
we found their spectral energy distributions consistent with pure photospheric emission. For IS Vir, we detect a marginal 
mid-IR excess which hints towards a dust disk. The orbital and rotational planes of IT Com appear to be coplanar, contrary 
O || to previous findings in the literature. We apply a multiple frequency analysis technique to determine photometric periods, 

and possibly changes of periods, ranging from days to decades. New rotational periods of 21.55±0.03 d, 65.1±0.3 d, and 
23.50±0.04 d were determined for V2253 Oph, IT Com, and IS Vir, respectively. Splitting of these periods led to tentative 
detections of differential surface rotations of 5P/P ~ 0.02 for V2253 Oph and 0.07 for IT Com. Using a time-frequency 
technique based on short-term Fourier transforms we present evidence of cyclic light variations of length ~ 10 yrs for 
V2253 Oph and 5-6 yrs for IS Vir. A single flip-flop event has been observed for IT Com of duration 2-3 yrs. Its exchange 
of the dominant active longitude had happened close to a time of periastron passage, suggesting some response of the 
magnetic activity from the orbital dynamics. The 21.55-d rotational modulation of V2253 Oph showed phase coherence 
also with the orbital period, which is 15 times longer than the rotational period, thus also indicating a tidal feedback with 
. the stellar magnetic activity. 
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1 Introduction the use of plate archives spanning decades or even a century. 

. As examples we mention that Phillips & Hartmann (1 19781 1 

During the last two decades numerous small and specialized determined first time long-term starspot activity of BY Dra 

surveys were carried out in order to identify new cool stars and CC Eri from Harvard plates from half a century between 

across the H-R diagram that show photospheric and/or chro- 1900-1950, and Frohlich et al. (120021 [2006) used the Son- 

mospheric signs of magnetic activity (e.g. Lockwood et al. neberg plate collection for studying the long-term behaviour 

119971 Baliunas et al. 119951 Henry et al. |1996l Strassmeier et of the solar-type single active star EK Dra and the K-giant 

al. 12000 , Henry et al. 120001 De Medeiros et al. 120021 Wright primary of the close binary HK Lac, respectively. Vogt et al. 

et al. 120041 White et al. 120071 Strassmeier et al. 120121 . The (120041 1 found several new variables on the Sonneberg plate 

next and still ongoing step was to obtain more accurate ba- collection of 34 years, most of them are early type stars 

sic astrophysical parameters for these stars and to conduct from late-B to A, showing slow long-term changes on the 

time-resolved photometric and spectroscopic monitoring in timescale of a few thousand days. Recently Balona (2012) 

order to determine rotational periods and orbital elements . found rotational modulation on A stars from Kepler data 

The seemingly most interesting objects were followed up which might originate from starspots as one possibility, and 

with dedicated campaigns with automatic telescopes at even if this would be verified than longer-term magnetic variabil- 

higher cadence in order to sample their starspot evolution ity could also be behind the slow changes shown by Vogt 

(e.g. Olah et al. [19971 Jarvinen et al. [20081 Berdyugina & et al.'s (120041 ) early type new variables. The DASCH (Dig- 

Henry [2007l Strassmeier et al. 2008. Savanov & Strassmeier ital Access to a Sky Century at Harvard) project recently 

120081 Korhonen et al. 120 101 Lehtinen et al. 120121) . resulted in the discovery of three K-giants' long-term vari- 

A valuable tool to extend the time-base of the datasets ability by Tang et al ([20101), which are very probably active 

for studying the long-term behavior of magnetic activity is stars, similar to those studied in this paper. 



* Corresponding author: e-mail: olah@konkoly.hu 



©WILEY H ONLINE LIBRARY 



©2013 WILEY-VCH Verlag GmbH & Co. KGaA. Weinhcim 



790 



Olah et al.: Long-term photometry of three active red giants in close binary systems 



Our main aim is to find evidence for a relation of pho- 
tospheric starspot activity with the evolutionary state of the 
star or the orbital properties of the binary system. Does the 
rearrangement of the internal structure at the base of the 
red-giant branch reflect onto changes of the dynamo prop- 
erties? Or, can tidal interaction dissipate orbital energy into 
magnetic energy? In this paper, we report on follow-up ob- 
servations of three chromospherically active binaries, each 
containing a K-giant star. 

The spectral types, effective temperatures, and radii 
(when known) of the primary stars in these binaries are 
similar to each other. However, their degree of rotational 
synchronism to the orbital revolution is markedly different. 
V2253 Oph is a short-period asynchronously-rotating star in 
a binary with an orbital period 15 times longer than the stel- 
lar rotation period. Its orbit is only slightly eccentric, if at 
all, because zero eccentricity is within 3er of the value given 
by Fekel et al. (2002). Thus, the period difference could not 
be due to pseudosynchronization. IT Com is also an asyn- 
chronously rotating star but its rotation period is longer than 
the orbital period by about 10%. Its inclination of the rota- 
tional axis is coplanar to the orbit, which itself is highly 
eccentric (Griffin 1988 1. The third system, IS Vir, is a syn- 
chronized giant in a marginally eccentric orbit (Fekel et al . 
120021 

In the following, we first summarize the relevant litera- 
ture for the three binary systems in Sect. |2] In Sect. [3] we 
describe our new observations and present the tools that will 
later be applied in Sect. [4] and Sect. [5] Sect. [6] summarizes 
and discusses our results. 

2 The program stars: V2253 Oph, IT Com, 
and IS Vir 

Table Q] summarizes some of the relevant basic parameters 
for the three systems for quick reference. 

The light variability of V2253 0ph (HD 152178) was 
discovered by Hooten & Hall (1990) who determined a 
rotational period of 22.35±0.05d from photometry. Their 
observations were initiated by the spectral classification of 
Houk d 19821 1 who noted Ca II H&K emission and found the 
star to be a single-lined (SB1) spectroscopic binary. She 
also derived a spectral type of G8/K0 (but erroneously of lu- 
minosity class V). The correct spectral type of the primary 
and the first orbital elements were published in the second 
CABS catalog (Strassmeier et al. 1993 1 showing a signifi- 
cantly longer orbital period of 314.5 d. The most recent or- 
bital elements are given in Fekel et al. (1999) and basically 
confirmed the original orbit. Strassmeier (11994) and Strass- 
meier et al. (U9941 determined absolute emission-line fluxes 
of the strong Ca II H&K emissions and ROSAT detected the 
X-ray flux of the system (Dempsey et al. 1 19931 >. 

IT Com (HD 118234) was first mentioned as a spectro- 
scopic binary by Griffin (1988) who published an SB1 or- 
bital solution with a high eccentricity (e=0.589). Together 
with unpublished photometry, he concluded that the primary 



should be an early K giant. Strassmeier (119941 1 and Strass- 
meier et al. ([1994) had this object on their target list and 
detected strong Ca II H&K emission lines while ROSAT de- 
tected the system in X-rays (Dempsey et al. 1993 ). The pho- 
tometric variability with a period of 64±1 days was discov- 
ered by Henry et al. (1995). This period implied that IT Com 
rotates slightly asynchronous with respect to the orbital pe- 
riod of P orb =59.054d. Henry et al. (119951 1 also noticed the 
strongly variable photometric amplitude between zero and 
0™2. Glebocki & Stawikowski (119971 ) stated that the orbit 
is not coplanar, i or b = 13° while i Iot = 34°; the authors 
found a general misalignment of the primaries of the asyn- 
chron systems to their orbits from observations and statisti- 
cal analysis. However, the error of the inclination determi- 
nations are between 5-9 degrees for systems like IT Com, 
which may question some of the individual results. Finally, 
IT Com was observed with the ELODIE spectrograph at 
OHP in the survey of Soubiran et al. (2003), who derived 
atmospheric parameters, distance, and space velocities to- 
gether with a large sample of mostly clump giant stars. 

The third target, IS Vir (HD 1 13816) is the best studied 
of the three targets in this paper. Fekel et al. (20021 summa- 
rized the various measurements of this object from X-rays 
to optical photometry and spectroscopy that was published 
until 2002. They obtained a new, well-determined orbital 
solution with synchronous rotation and found a very low or- 
bital and rotational inclination of this SB1 binary. From 12 
years of photometry, cyclic light variability of the order of 
7-8 yrs was discovered. Detailed photospheric abundances 
of IS Vir were published by Katz et al. (120031 ) and by Morel 
et al. ( 120031 1. Both analysis found near-solar abundance with 
perhaps a slight overabundance of Na, Mg, Al and Ca with 
respect to Fe, by 0.1-0.2 dex; but the derived Fe/H ratio is 
quite different (0.04-0.09 and -0.11 ± 0.09, respectively. 

3 Observations and data-analysis methods 
3.1 New APT observations 

We present altogether 3,191 new photometric Vic data 
points. All observations were taken with the 0.75m 
Vienna-Potsdam Automatic Photoelectric Telescope (APT) 
Amadeus at Fairborn Observatory in southern Arizona 
(Strassmeier et al. [T997T >. For V2253 0ph, 1,883 measures 
were available from between June 1997 and June 201 1. For 
IT Com, 698 measures from December 1999 through April 
2010 and for IS Vir 610 measures from December 1999 
through January 201 1 could be used. These data are plotted 
in Fig. Q] for reference. Observations of V2253 Oph from 
the first half of 1997 were published earlier (Strassmeier 
et al. 1999) but are added to the present dataset. An inde- 
pendent continuous 12-yr photometric dataset for IS Vir was 
published and analyzed by Fekel et al. (2002) who started 
observing 10 years prior to our APT campaign with two 
seasons overlapping. We combine these two data sets and 
use it for parts of the present analysis. Our Amadeus APT 
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Fig. 1 Overview of the Vienna-Potsdam APT light curves of the three binaries in this paper; a. V2253 Oph, b. IT Com 
and c. IS Vir. The APT took at least one Vic data point every clear night for the past 15 years. In the figure, target stars 
are arranged in columns. V magnitude is plotted in the top row, Ic magnitudes in the middle row, and V — Ic color in the 
bottom row. 



photometry had larger than usual scatter in the years be- 
tween 2007 and 2009 due to target acquisition problems 
(see Strassmeier et al. 120121 for the details). However, the 
typical external precision of Amadeus has been around A— 
6 mmag depending on the brightness of the target. Photom- 
etry from literature sources are sparse for the other two sys- 
tems, V2253 0ph and IT Com, and are not taken into ac- 
count because we would scrutinize the homogeneity of the 
APT data. Note that the IS Vir data of Fekel et al. J20021 
were taken with the neighboring TSU 0.4m T3 APT at Fair- 
born Observatory (but with different comparison and check 
stars) and its basic data reduction coefficients are therefore 
practically identical and the two data sets straightforward to 
combine. 

All data were taken differentially with respect 
to a nearby comparison star and a check star. A 
nearby bright star is used for (automatic) navigation. 
The comparison/check stars were HD152501/HD151 179, 
HD118670/HD 118244, and HD1 13448/HD1 12805 for 
V2253 Oph, IT Com and IS Vir, respectively. A single mea- 
surement is the mean of three readings between the variable 
and the comparison star. Sky readings through the same 30" 
aperture as for all other readings are taken usually between 
the comparison and the variable star measure to ensure time- 
liness. The detailed data reduction and Johnson-Cousins 
transformations for Amadeus were described in Granzer et 
al. (1200 IK We refer the reader to this and subsequent papers 
as well as references therein. 



3.2 Determining photometric periods 

Throughout this paper, we use the program package 
MuFrAn (Multiple Frequency Analysis; Kollath 119901 ) to 
quantify rotational periods and its errors, which are deter- 
mined increasing the residual scatter of the least-squares so- 
lutions by about 10% of the precision of the data. Another 
package TiFrAn (Time Frequency Analysis; Kollath & Olah 
(2009) allows a time-frequency analysis with several meth- 
ods if the datasets are long enough in order to check for 
possible cycles of length of a few years. In this paper we 
apply a Short-Term Fourier Transform (STFT), numerical 
details of the procedure are given in Kollath & Olah (2009), 
including the input of data with variable (rotational) modu- 
lation and in the presence of seasonal gaps. This algorithm 
was already used for the detection of cyclic behavior of sev- 
eral other active stars (e.g. Olah et al. 2009 1. 



3.3 Fitting spectral energy distributions (SEDs) 

We determine effective temperatures by fitting ATLAS -9 
models (Castelli & Kurucz 2003) to optical and infrared 
fluxes. Optical fluxes at Johnson BV and Cousins Ic wave- 
lengths were derived from our light curves for the times of 
maximum stellar brightness (presumably the least spotted 
state). In order to take into account the large variability of 
stars, we assigned an uncertainty of 0. 1 mag to each photo- 
metric value. 
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Table 1 The pro, 


gram stars. 






Parameter 


V2253 Oph 


IT Com 


IS Vir 


HD 


152178 


118234 


113816 




8 5 

O.J 


7 4 


O. J 


M-K class 


KOIII 


KOIII 


K2III 


1 eft y^) 


4900+ 1 00 


4700 


4700±77 


iu & y 


(2 5^ 


3.1 


2 45+0 21 


n Qin 7 ^lcm s ^ 


28.8±1 


5.3 


5.9+1 


1 L oil! £ ^1\.(V) J 


12.7 


7 


2.7 


1 \ ) 


63 


13 


13+4 


[Fe/H] 




—0.19 


—0 11 + 09 


Orbit 


SB1 


SB1 


SB1 


P orb (days) 


314.47±0.18 


59.054±0.004 


23.655+0.001 


e 


0.024±0.009 


0.589±0.010 


0.022+0.008 


K(km s _1 ) 


14.32±0.15 


9.66±0.15 


6.598+0.054 


f(m) 


0.096+0.003 


0.0029±0.0002 


0.00070+0.00002 


ai sini (10 6 km) 


61.92±0.64 


6.34±0.12 


2.146+0.017 


d (pc) 


360±105 


185±25 


250+50 



Source of orbital elements: Fekel et al. (IT9991 for V2253 Oph, Griffin ( 1988 1 for IT Com and Fekel et al. d2002l for IS Vir. The distance 
is from the revised Hipparcos parallax (van Leeuven 20071. The other parameters are for orientation only and were taken from the 
CABS-III catalog (Eker et al. 2008) and references therein. 



Near-infrared JHK S data were taken from the Two Mi- 
cron All-Sky Survey (2MASS; Skrutskie et al. 120061 These 
data were supplemented by photometry from the Wide-field 
Infrared Survey Explorer (WISE) satellite. We used the Wl 
band (centered on 3.4 /im) photometry from the All-Sky 
Database (Wright et al. 120101 1 but rejected the M^2-band 
(4.6 /^m) measurements because these were saturated for all 
of our targets. 

We fixed both the surface gravity (log g) and the metal- 
licity ([Fe/H]) during the fitting procedure. For IT Com and 
for IS Vir these values were set to the literature data cited 
in Table [T] (Morel et al. |20031 l. We then selected model 
atmospheres from the ATLAS-9 grid that have the clos- 
est metallicity and log g value to the quoted ones. In the 
case of V2253 Oph no literature data for gravity and metal- 
licity were found. Therefore, we adopted solar metallic- 
ity and estimated the surface gravity from its spectral type 
and luminosity class (see Table All of our targets have 
good quality photometry in the mid-infrared Wi (12 /im) 
and WA (22 /im) bands with signal-to-noise ratio >20:1. 
By utilizing the W3 and WA bands, we are also able to 
search for dust emission around our targets. None of our 
sources are affected by any known artifacts. As a first step, 
the predicted photospheric flux densities were derived for 
the WISE wavelengths using the best-fit ATLAS-9 mod- 
els. Then the significance level of any infrared excess was 
calculated for both mid-infrared WISE bands by xa = 
(-Fwise — -Fpred)/cA, where Fwise is the flux in a specific 
WISE band and F prc d is the corresponding predicted flux 
value from ATLAS-9. a\ is the quadratic sum of the uncer- 
tainty of the measured flux density, the absolute calibration 
uncertainty (taken from the Explanatory Supplement of the 
WISE All-Sky Data Release Products), and the uncertainty 
of the predicted flux density in the specific band. Following 



the general approach for clear excess detection, we require 
Xx > 3.0 (e.g. Kennedy & Wvatt l20T2l 



3.4 Estimating stellar ages 

Ages are estimated from a comparison of the location in the 
H-R diagram with evolutionary tracks from Pietrinferni et 
al. (|2006). Stellar luminosity is based on the Hipparcos par- 
allax (van Leeuwen 2007) and the brightest magnitudes of 
our y-band light curves. Bolometric corrections were com- 
puted based on T c h values using the formula outlined in Tor- 
res (|2010), and we adopted an absolute bolometric magni- 
tude of 4.73 for the Sun. Luminosity errors are simply prop- 
agated parallax and magnitude errors. The temperatures and 
its errors are those from the SED fits. 

There is an indirect dynamical method to estimate the 
age of a binary system, or at least to give some constraints 
for it. The method is based on a comparison of the present- 
day orbital configuration with the theoretical tidal circular- 
ization time scales. At this point we refer to Verbunt & Phin- 
ney d 19951 ) who investigated the circularization in binaries 
where the primary component is an evolved, red giant star. 
From the equilibrium-tide theory of Zahn (11977119 89) these 
authors derived a practical formula (their Eq. [7]) which 
gives the rate of circularization for every evolutionary stage 
of an evolved primary as a function of the mass-ratio and 
the orbital period. We note that this formalism depends only 
weakly on the (uncertain) mass-ratio but strongly on the ac- 
cepted evolutionary model, and thus on the mass of the pri- 
mary. 
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Table 2 Revised stellar parameters. 



Wavelength [urn] 



Fig. 2 Spectral energy distributions of V2253 0ph (top), 
IT Com (middle) and IS Vir (bottom). The dots are the ob- 
served fluxes from our APT data, from 2MASS, and from 
WISE. The lines are the best fit ATLAS -9 models. Effective 
temperatures and interstellar extinction values are given in 
each panel. 

4 Masses and ages 

The spectral energy distributions (SEDs) of V2253 Oph, 
IT Com and IS Vir are plotted in Fig. [2] and the result- 
ing effective temperatures and extinction values are listed 
in Table[2] Detailed H-R diagrams are given in Fig. [3] 



4.1 V2253 Oph 

For V2253 Oph, we found the near-IR and mid-IR fluxes 
to be consistent with pure photospheric emission of T e g = 
5050±200K. This value is close to the recent determination 
of Bailer- Jones d20TTl of 4944 and 4814^25^, using 
BVJHK photometry and parallax with two different meth- 
ods. Our quite large E(B-V)=0.24±0.08 is also supported 
by Bailer-Jones' d20TTb values of 0.18t^ 2 and 0.14±g;^ 
from the two methods, respectively, Errors mark the 90% 
upper- and lower bounds for both parameters. 



Parameter 


V2253 Oph 


IT Com 


IS Vir 


Vmax 


8.61 


7.45 


8.30 


Vmiii 


8.87 


7.86 


8.70 


(V - Ic) 


1.25-1.35 


1.12-1.24 


1.17-1.28 


P ro t (days) 


21.55±0.03 


65.1±0.3 


23.50±0.04 


Cycle? 


10 yr 


flip-flop 


5-6 yr, 2 yr 


T cff (K) 


5050±200 


4750±130 


4600±90 


E{B - V) 


0.24±0.08 


0.05±0.05 


0.07±0.04 


logl/i (L©) 


1.97±0.37 


1.68±0.15 


1.66±0.21 


Ri (R©) 


12.7±5.4 


10.3±1.8 


10.7±2.6 


Mi (M r 


«3 


«1.4 


wl.3 


M 2 (M ) 


1.35: 


0.29: 


0.55: 


in 


^65 


40±12 


15±5 


oi (10 6 km) 


68.3 


9.38 


8.87 


Constant values used when calculating the revised parameters 




2.5 


2.5 


2.5 


Fe/H c 


0.0 


-0.2 


0.0 


distance (pc) 6 


360±105 


185±25 


250±50 



11 Possible mass ranges are given in the text. For IT Com and 
IS Vir Fe/H is from the literature, for IS Vir we suppose solar 
abundance because of the contradicting (but close to the solar) 
literature values. The logg of IT Com of Soubiran et al (20031 
has no published error and is considered uncertain. c From the 
corrected Hipparcos parallax. 



The position in the H-R diagram in Fig. [3^ shows that 
V2253 Oph is quite massive, about 3 M Q with a mass range 
between 2.4-3.8 M Q . Hence, the system cannot be very old, 
maybe 1 Gyr or even less according to the models of Pietrin- 
ferni et al. d2006l . 

From the formula of Verbunt & Phinney (119951 we find 
that in case of the «314-d orbital period, we could expect 
a circularized orbit only as late as during the AGB phase 
of stellar evolution. Apparently, the star is not in a such 
late phase of its evolution. A possible solution is offered 
by Verbunt & Phinney (1995). They suggest, that for some 
circularized very wide binaries, a possibility to explain the 
circular orbit with such a massive primary is, that the un- 
seen secondary could be a white dwarf. It is thought that the 
white dwarf forces the circularization during the primary's 
red giant phase either due to the stronger tidal friction dur- 
ing the epoch when the primary radius becomes very large 
or/and a more rapid mass-transfer in case of a Roche-lobe 
filling configuration. 

The orbit of the V2253 Oph system is fairly wide, so 
the suggestion of Verbunt & Phinney ( 1995) could be con- 
sidered. The minimum radius of the primary of 12.7 R Q to- 
gether with the rotational period of 21.55 d and a measured 
vsrai of 28.8 kms" 1 (Fekel et al. [T9991 and 24 kms" 1 
(Strassmeier et al. 1994) yield moderately high inclinations 
of the rotational axis of i ~ 70° and 54°, respectively for 
these v sin i values. Assuming a coplanar orbit, this suggests 
a mass ratio between q^0.43-0.51 resulting in a secondary 
with a mass of about 1.3 M Q orbiting the 3 M Q primary. 
On the other hand, if the primary mass is just 2.5 M©, then 
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Table 3 Newly determined galactic coordinates and ve- 
locities 



par. (unit) 


V2253 Oph 


IT Com 


IS Vir 


v rad (km/s)° 


-36.75±0.13 


-19.13±0.08 


21.25±0.04 


U(km/s) 


-34.9±0.7 


-63.8±9.0 


12.4±1.5 


V(km/s) 


-13.8±7.1 


-88.6±13.2 


-26.8±4.2 


W(km/s) 


-19.9±5.7 


-5.3±2.1 


7.2±2.5 


X(pc) 


352.1 


37.8 


85.1 


Y(pc) 


-29.9 


1.9 


-101.5 


Z(pc) 


67.5 


180.6 


210.5 


b (deg.) 6 


10.81 


78.18 


57.82 



"Gamma velocities form the spectroscopic orbits. 'Galactic lati- 
tude. In the calculation of the Galactic space velocity we used a 
right-handed coordinate system (U is positive toward the Galactic 
center, V is positive in the direction of Galactic rotation, and W 
is positive toward the north Galactic pole) and followed the gen- 
eral recipe described in the Hipparcos and Tycho catalogs (ESA 
1997). X, Y, Z are the physical space coordinates, centered on 
the Sun, in the same directions as U, V, W.For the calculation 
of galactic space motion components the coordinate, proper mo- 
tion and trigonometric distance information were taken from the 
Hipparcos catalogue (van Leeuwen 2007 ). 

the secondary mass would be slightly less than 1.2 M Q . 
More about the mass of the secondary is in Sect. [B] in the 
appendix. 

The existence of such a high mass white dwarf has been 
proven by Kalirai et al. (2005), who studied white dwarfs 
in the young (650 Myr) cluster NGC 2099. For WD 24 
a present-day mass of 1.11 M Q was claimed and an ini- 
tial mass of 4.43 Mg computed, although with consider- 
able errors. Thus, a scenario of V2253 Oph with a red-giant 
primary of 2.4 M Q and a white-dwarf secondary of about 

1.1 M Q , which should have been originally over 4 M Q , is 
not impossible, taking into account the young age of the sys- 
tem, which could have been still long enough for the present 
secondary to evolve to a high-mass white dwarf. 

4.2 IT Com 

For IT Com, we found the measured mid-infrared fluxes 
to be consistent with pure photospheric emission of T e g = 
4750±130K. Bailer-Jones' d20TH ) results of 4821+^ and 
4735I213K with two different methods are practically the 
same as our result, although with quite high uncertainty. The 
near zero E(B-V)=0.05±0.05 we found is lower than that of 
Bailer- Jones's dlPTTT l of 0.14t°;^ and 0.11 ±0.09 from the 
two methods, but still are well within the respective errors. 

The position in the H-R diagram suggests a compara- 
ble low mass of 1.4 M© (Fig. |3})) with a full range of 1.2- 
1.6 M Q . The Pietrinfemi et al. (2006 1 tracks indicate an age 
of ss 1.6 Gyrs. Large dispersion in the galactic space veloc- 
ities of IT Com (Table [3) implies that it may not be a thin 
disc star. Indeed, using the kinematic criteria proposed by 
Reddy et al. (2006 ) we found that IT Com probably belongs 
to the thick disc population. 

© 2013 WILEY- VCH Verlag GmbH & Co. KGaA, Wcinheim 
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Fig. 3 Positions of the three primary stars in the H-R 
diagram (dots with error bars). From top to bottom, a. 
V2253 Oph, b. IT Com, and c. IS Vir. The vertical axis plots 
logarithmic luminosity in solar units and the horizontal axis 
is logarithmic effective temperature in Kelvin. The lines are 
tracks for different masses as indicated in the inserts and 
were computed for solar metallicity for V2253 Oph and 
IS Vir, but interpolated to [Fe/H]=-0.2 for IT Com. 



The orbit has a high eccentricity which, together with 
the dynamic age based on Verbunt & Phinney (1995 ), sug- 
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gests that IT Com is in some stage of H-shell burning and 
not yet on the AGB because then the orbit should have al- 
ready been circularized. 

The calculated radius of 10.3 Rq and wsini of 
5.5±1 kms -1 (Henry et al. 1995 ) constrain the inclination 
of the rotational axis to between 34-54°, in good agreement 
with the most likely coplanar orbital inclination of about 
40°, derived in Sect.|A]in the appendix. The mass of the sec- 
ondary star is about 0.55M©, see Sect. |B]in the appendix. 



4.3 ISVir 

For IS Vir, we derived T oS = 4600±90K, very similar to 
the recent determination of Bailer- Jones (201 1) of 4521 
and 4574±^?K. Our E(B-V)=0.07±0.04 is practically the 
same as that of Bailer-Jones' d20TTt of 0.03±g;g£ and 
0.06 1° 06 from the two methods. 

The SED fits for the mid-IR bands achieved XW3 = 1-5 
and xwa = 2.7, implying a marginal excess detection in the 
22-/im W4 band. Recently, a comparable amount of excess 
emission at WA was announced for HD 139357, a planet 
hosting K4 III giant star (Morales et al. 120121 1. Generally, 
Jura d2004b explained excess IRAS 25 /im emission by dust 
release during ice sublimation of Kuiper belt-like objects as 
the star evolves from the main sequence onto the red giant 
branch. Further infrared observations are needed to confirm 
the excess and to clarify whether it is related to circumstel- 
lar dust or can be explained differently, e.g. by the contri- 
bution of a background galaxy or by an apparent positional 
coincidence with an interstellar dust cloud. 

The position of IS Vir in the H-R diagram is compatible 
with a mass of 1.3 M (Fig. [3};), i.e. of similar mass than 
IT Com. The mass range is 1.1-1.5 and the age of the system 
is certainly beyond 5 Gyrs from the HRD. 

Fekel et al. ( 2002 ) derived a similar mass of 1 .5 M e and 
claims the system is not young since the star is about 250 pc 
above the galactic plane. The tidal circularization time scale 
does not provide any further constraints for its age because 
the orbital-period limit for circularization is longer or com- 
parable (28 d) to today's orbital period (23.6 d). 

IS Vir has a low orbital inclination of about 13° accord- 
ing to (Fekel et al. 2002). From its newly calculated radius 
of 10.7 R Q , its rotational period of 23.5 d, and the projected 
rotational velocities found in the literature. The inclination 
lies between 10-20°; the range of inclination comes from 
the uncertainties of the radius and the different velocity data 
collected in the CABS-III catalog (Eker et al. 120081 . The 
star is seen nearly pole-on. As a consequence of the low 
inclination the amplitude of the rotational modulation is al- 
ways small and less than 0^1 in V, as can be seen in Fig. QJ. 
The mass of the secondary component must be w0.55 M , 
see Sect.|B]in the appendix, for more. 



V2253 0ph P rot= 2 , L55d Fourier spectrum 




0.04 0.045 0.05 




-0.005 0.005 



Frequency [c/d] 

Fig. 4 MuFrAn periodogram of V2253 Oph from the V- 
band data until the end of 2005. Two significant periods are 
marked in the upper panel, the spectral window is plotted 
below. 

5 Starspot activity 
5.1 V2253 Oph 

The F-band amplitude due to rotational modulation var- 
ied between 0l n 15 in 1997 at a time of faint average light 
and 0' n 10 in around 2002 at a time of bright average light 
(Fig.QJ). This suggests a relation between the overall spot- 
tedness and the size of individual spots or spot groups sim- 
ilar to what has been seen on other over-active stars. How- 
ever, the shape of the light curve of V2253 0ph changes 
continuously and smoothly between subsequent rotations. 
This is indicative of rapid spot rearrangements typically 
observed for the more active stars. Not surprisingly, the 
MuFrAn periodogram reveals two closely-space periods, 
but only until the end of 2005, see Fig. [4] The dominant 
21.55±0.03-d period with an average 0™06 full V am- 
plitude is seconded by a 22.05±0.04-d period with a full 
0™025 amplitude. We tentatively interpret this due to dif- 
ferential surface rotation and estimate a relative shear factor 
SP/P of «0.02, i.e. a lap time of the equator with respect 
to the poles of «314d. 

The rotation of the spotted primary is strongly asyn- 
chronous with respect to its 314-d orbital period. Note that 
the eccentricity is very low and does not exceed its 3 a un- 
certainty, and can not account for the period difference. Yet 
the light variation of two subsets of our APT data seem to 
produce a fairly homogeneous light curve when plotted ver- 
sus the orbital period (Pig. [5J. Its amplitudes are marginally 
larger than the "scatter" but show extrema (maximum and 
minimum) always near periastron. The very first seasonal 
data set from 1997, which shows the highest rotational- 
modulation amplitude of all of our data, appears to have a 
minimum at orbital phase near 0.5, half between the formal 
times of periastron and apastron (note that one observing 
season usually covers 140 nights, i.e. one half of the orbital 
period). Possibly, the 1997 season belonged to yet another 
orbit-induced variability pattern. Due to the small eccentric - 
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Fig. 6 Time-frequency diagrams for V2253 Oph. a. from the V magnitudes and, b., from the V — Ic color index. The 
top panels show the observations as red (grey) dots, blue (black) dots show the data prewhitened with the main rotational 
period and its harmonics, red and blue (grey and black) lines show spline interpolations to the respective data. Notice the 
strongest signals are at around 10 yrs from V and «2.5 yrs from V — Ic- 



ity it is not very likely that this variation originates from 
magnetic interaction of the two binary components at peri- 
astron but nevertheless there remains some evidence for a 
relation of the spot activity on the primary with the orbit of 
the system. 

The time-frequency diagram from 15 years of photom- 
etry of V2253 Oph is plotted in Fig. [6] The most prominent 
long-term variation is even seen by eye and stretches over 
10 yrs with a peak-to-peak V amplitude of (f. n 1 . Despite that 
the time-frequency diagram gives a strong signal around 10 
years, we can not be sure of its cyclic nature due to the rela- 
tively short time coverage. Half of this period, a variability 
of about 4.5 years, also appears in the time-frequency dia- 
gram. The similarly constructed plot from the V — Ic color 
index does not show any traces of the dominant 10-yr vari- 
ation but shows a weak signal around 4.5 years. However, 
the largest amplitude in the V — Ic color index variation 
is reconstructed on the timescale of 2.5 years. All these pe- 
riodicities seem to be related to each other because within 
their uncertainties of about 0.5-1 yr they appear to be har- 
monics. Note that the color scales in the two plots in Fig. [6] 
are the same, which allows to directly compare the ampli- 
tude of the signals. We also note that despite the scatter of 
the measurements was larger between 2007-2009 due to an 
instrumental problem, the comparison-check data proofed 
that the mean light level was not affected, so we conclude 
that the large long-term variation is real. 



Unfortunately, we are left with an inconclusive situa- 
tion. It is tempting to interpret the long-term V^-brightness 
variation due to a dynamo-induced change of the over- 
all spottedness, like for the solar 11 -year cycle, and the 
long-term V — Ic index variation due to an average spot- 
temperature change. However, the two tracers should vary 
with the same "cycle" period or, alternatively, are not re- 
lated with each other. A bit frustratingly we conclude that 
after 15 years of photometry we need another 15 years of 
photometry. 

5.2 IT Com 

The best-fit average photometric period for IT Com from 
the 11-yr U-band data is 65.1±0.3d (Fig.[7|l. However, the 
photometry shows significant light-curve changes from one 
stellar rotation to the next (note that only two stellar ro- 
tations can be covered during one observing season) that 
may affect its period determination. An amplitude increase 
of the rotational modulation by nearly a factor of three hap- 
pened between 2005 and 2006. In 2002, a sudden phase shift 
by «0?4 was witnessed during the observing season which 
then reverted back to the original phase in early 2005. This 
could be termed a classical flip-flop phenomenon (e.g. Elst- 
ner & Korhonen|2005] Olah et al. 120061 . It has stalled near 
the periastron passage of the two binary components at JD 
2,452,406.734, marked in Fig.[8]top panel where the obser- 
vations before, during, and after the phase flips are plotted 
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Fig. 5 Evidence for a relation of the spot activity on the 
primary of V2253 Oph and the binary orbit. The top panel 
shows again the long-term AV magnitudes plotted versus 
Julian date while the bottom panel shows the same data plot- 
ted versus orbital phase (P rb = 3 14. 17 d). Symbols on both 
panels refer to the same data points. 
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Fig. 8 Upper panel: AV magnitudes of IT Com versus 
JD for 2000-2010. The plot symbols and colors mark the 
data embraced by the flip-flop. The first phase flip started in 
early 2002 and flopped back by early 2005. Lower panel: 
The same data as above but phased with the rotational pe- 
riod of 65.1 d. The plot symbols and colors are the same 
as in the top panel and indicate the phase shift of the data 
within the flip flop. Note that the lower panel normalized the 
differential magnitudes to the mean light level of the three 
subsets. 



Fourier spectrum 




Frequency [c/cTj 



Fig. 7 MuFrAn periodogram of IT Com from the V"-band 
data. Three significant periods are marked in the upper 
panel, the spectral window is plotted below. 



with different symbols (and color). The back-flip in 2005 
was actually not directly seen because it happened outside 
the annual observing window but the phase coherence there- 
after was the same as before the flip and constant until the 
end of the APT observations in 2010. The time duration 
from flip to flop was between 2.0-2.7 years, i.e., it ended 
between the last observation of 2004 and the first one in 
2005. So far, the phenomenon did not repeat and we could 
not speak of a flip-flop "cycle" but state that flip flops may 
occur aperiodically and its length is not related to its reoc- 
currence. 

Systematic seasonal changes are also obvious for the av- 
erage U-light level in Fig. [Tjb (replotted on a larger scale 
in Fig. [8] top panel) and amount to almost 0". n 15. Unfortu- 
nately, the 1 1 years of photometry is insufficient to search 
for cyclic changes; a simple look at the data suggests a 
timescale of a cycle (if any) longer than the length of the 
dataset itself. The long-term V — Ic light curve traces the 
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Fig. 9 Phased light and color index curves of IT Com with 
the 65.1 days period, phase zero is at periastron. The plot 
symbols are the same as in Fig. [8] V — Ic is redder when 
the star is fainter. 

average V^-light level in the sense that redder V — Ic oc- 
curs when the system is brighter. The V — Ic variation is 
in phase with the color variability, at the same time the am- 
plitude of the rotational color modulation does not relate to 
the amplitude of the V (or Ic) modulation, which indicates 
a constant spot temperature, see Fig. [9] 

Apart from the dominant rotational signal of 65.1 d 
(with (J} 1 1 1 full V amplitude), two more significant peri- 
ods show up, see Fig. [7] One is slightly longer, 66.6 d with 
an amplitude of 0Y 1 072, and one slightly shorter, 63.9 d with 
an amplitude of I . n 068. Both periods have an approximate 
error of ±0.3 d and are highly significant. 

Between 2000-2002 the 65.1-d period is dominant (red 
dots in Fig.[8]l, whereas during the next two years the 66.6-d 
period gave the best fit (blue crosses in Fig. |8). From 2005 
onwards again the 65. 1-d period was dominant and resulted 
in the least scattered light curve. This indicates that the indi- 
vidual periods stem not from longitudinal migration due to 
some sort of a differentially rotating surface but are related 



with the emergence, or decay, of an active longitude as part 
of the flip-flop phenomenon. If we nevertheless assume a 
solar-like differential rotation law, adopt a latitude range for 
the spots of 50°, and assume the shortest period belongs to 
the equator and the longest to a latitude of 50°, we may es- 
timate a lower limit for the relative differential rotation of 
5P/P « 0.07. Therefore, the lap time of the equator with 
respect to the high latitudes is «90 d. 

5.3 IS Vir 
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Fourier spectrum 
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Fig. 10 MuFrAn periodogram of IS Vir from the V^-band 
data. One period is marked in the upper panel, the numer- 
ous peaks around it indicate other, slightly different periods. 
Note the low amplitude of the variations, see the text for 
more. The spectral window is plotted below. 

Fekel et al. (120021 published 12 seasons of BV pho- 
tometry from the TSU T3 APT of which the last two sea- 
sons overlap with our first two of the 1 1 seasons of Vienna- 
Potsdam APT data. Based on the comparison of the 76 com- 
mon nights during the two overlapping seasons, we com- 
bined the two datasets with a shift of 2.07 mag. in V (due 
to two different comparison stars) to a 21 -year time series 
for further analysis. After pre- whitening with the long-term 
high-amplitude variations with a 6250 days period and its 
four harmonics (which is a purely mathematical tool to get 
rid of the long-term variation) we find a mean photometric 
period of 23.50±0.04d from Johnson V that we interpret to 
be the stellar rotation period. It is very close but not equal 
to the 23.655±0.001 d orbital period (Fig[l0). 

During the first three and the last three seasons when the 
star was the brightest, the analysis does not give clear peri- 
ods. The data at medium brightness, i.e., without the first 
and last three seasons and the season of the minimum, re- 
sults in two closely-spaced periods of 23.64 and 23.50 days 
with an error of 0.3 days and a full Fourier amplitude of 
0.024 mag. Finally, the year of minimum light shows a pho- 
tometric period of 22.87±0.80 days with a full Fourier am- 
plitude of 0.028 mag. Considering the low inclination of the 
star (w 15°) this result suggests, that at maximum light 
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Fig. 11 Phased light curves of IS Vir with the 23.50 days 
rotational period. Red and blue (dots and crosses) symbols 
mark our observations and that of Fekel et al. (2002). Note 
the very low or sometimes zero amplitude of the rotational 
modulation and the long-term brighntess change of about 
0.4 magnitudes. 
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Fig. 12 Time-frequency diagram for IS Vir. The top panel 
shows the combined data from the TSU T3 APT and the 
Vienna-Potsdam Amadeus APT. Observations are plotted as 
red (grey) dots, blue (black) dots show the data prewhitened 
with the main rotational period and its harmonics, red and 
blue (grey and black) lines show spline interpolations to the 
respective data. 



only the polar region has some spots which cause no or 



very small rotational modulation, while during other times 
spots may rotate in and out of view at lower latitudes, with 
shorter period and higher amplitude at minimum. However, 
according to the errors of the periods and the very low am- 
plitudes we consider this only as a possible scenario. The 
amplitude change is also seen in V — Ic (Fig. which 
suggests that at least part of the long-term variations is due 
to a temperature change rather than an increase or decrease 
of the overall spottedness. Fig. QT| depicts the very small 
amplitude rotational modulation and the large change in the 
overall light level of the star during 21 years of observations 
(blue crosses: Fekel et al. (120021 1. red dots: this paper). 

The full APT time series in Fig. Q~2] indicates a long- 
term variation with a duration of 5-6 yrs which is slowly de- 
creasing with time, and varying amplitudes of up to 0™3 in 
V. A weakly defined timescale of about 2 yrs also shows 
a decrease from 2.3 to 2.0 yrs. These cycle lengths are 
not harmonics of each other, therefore their parallel change 
suggests a common origin. The same feature, i.e., paral- 
lel change of two cycles have already been observed in 
V410 Tau, EI Eri and UX Ari by Olah et al. (|2009l >. 

6 Summary and discussion 

We have investigated three active (non-eclipsing) SB1 bina- 
ries with early K giants as primaries. The estimated ages of 
the three systems are naturally rather uncertain, V2253 Oph 
is about 1 Gyr old or less, while IT Com and IS Vir seem 
to be older, over 2 Gyrs, likely a few Gyrs. Without more 
precise distances these estimates can not be significantly 
improved at the moment. More realistic ages were found 
for two other active binary systems with giant components; 
BE Psc («2.7Gyrs, Strassmeier et al. 120081 and CZ CVn 
(6-7 Gyrs, Strassmeier et al. 1201 1 ). BE Psc is an eclipsing 
binary and therefore the stellar parameters are known much 
better than in case of SB 1 systems, and CZ CVn is less than 
100 pc away with just a small extinction, thus its luminosity 
is better constrained than that of the systems in the present 
paper. 

V2253 Oph has the highest mass («3 M ) of the three 
active binaries and is also the youngest. Its orbit is nearly 
circular but the rotation is very fast and strongly asyn- 
chronous, the rotational period is almost 15 times shorter 
than the orbital one. According to Verbunt & Phinney 
(1995) this unusual scenario suggests either a white dwarf 
companion or an initially low eccentricity of the orbit. Apart 
from the rotational modulation the light variation shows 
a systematic change for several years during the 15 times 
longer orbital period with extrema (maximum or minimum) 
around periastron. The two photometric periods found in the 
data until 2005 are probably due to differential surface ro- 
tation which indicates a relative shear factor of SP/P of 
«0.02, and a lap time of the equator with respect to the 
poles of «314d. That this value is nearly identical to the 
orbital period may be coincidental, we have no explanation 
for it. The activity has a cycle with a characteristic length of 
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about 5 years and a.V — Ic color index, thus temperature, 
variation of about 2.5 years (i.e., half cycle), plus a long- 
term variation comparable to the length of the dataset. In the 
literature only two, barely studied similar systems with gi- 
ant primaries, low (zero) eccentricity and much shorter rota- 
tional than orbital periods, are found: DW Leo and VX Pyx, 
however, in both cases the rotational modulation is found 
be marginal thus their periods are uncertain (see CABS-III: 
Eker et al. |2008l for the details). 

IT Com and IS Vir have similar masses (1.4 and 
1.3 M Q ), and their orbits are of similar size (a! is 
9.38 xlO 6 km and 8.87 xlO 6 km, respectively), however 
the eccentricities are different: IT Com has a very eccen- 
tric orbit while IS Vir has nearly circular. Both systems are 
older than V2253 Oph. IT Com rotates 9% slower than its 
orbit and has a very low mass secondary (about 0.3Mq) 
so the tidal interaction is not very effective. On the con- 
trary, IS Vir which has a higher mass secondary (about 
0.55M Q ) thus more effective tidal interaction, is synchro- 
nized, and the shorter orbital period is below the period 
limit of circularization. The primaries have strong magnetic 
activity and the magnetic braking may slow down the ro- 
tation, below even the orbital period in case of IT Com, 
where the low mass secondary cannot effectively interact. 
Strassmeier et al. (120121 studied V383 Ser (HD 142680), an 
SB2 system, the only one which shows some similarity with 
IT Com: it has a considerable eccentricity (e=0.3158) and 
rotates super-synchronously (P or b=24.534 days and P ro t= 
33.4 days), but no long-term photometry is available for fur- 
ther studies. 

We directly observed an abrupt exchange of the deeper 
and shallower minima (stronger and weaker active longi- 
tude) of IT Com in 2002. The new configuration lasted for 
about two years, after that the previous scenario returned. 
During the 1 1 years of observations three photometric pe- 
riods were always present strongly suggesting a differential 
rotation with a minimum a ~ 0.07. From a combined, 22 
years long dataset of IS Vir we found evidence for a slowly 
decreasing cycle period between about 6-5 years. 
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A On the inclinations of IT Com 
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Fig. Al Possible orbital inclinations of IT Com with dif- 
ferent mass ratios. Red solid line and blue dotted line mark 
primary masses of 1.4 and 2.5 solar masses. 



The orbital elements of IT Com were derived by Griffin 
(1988), who found a very eccentric orbit (e=0.589) and a small 
mass function (/(m)=0.0029). Griffin suggested that the primary 
was a K0 or a Kl giant and the system very similar to 73 Leo 
consisting of a K giant primary and an Fl secondary. Based on 
this information, Glebocki & Stawikowski ( 1997 1 adopted a mass 
of 2.5 M© for the K giant primary, a mass ratio of 0.64, and from 
this got the low orbital inclination of 13°. This calculation used 
the fact that when the mass function is extremely low then the 
orbital inclination is not very sensitive on the primary mass and 
the mass function. The relation between the orbital inclination 



sin* = Kl Vl^i {/ (1 + g)3 £ , (Al) 
q V Gmi 

where K\ is the amplitude of the radial velocity curve, e is the 
eccentricity, mi is the mass of the primary, q = 7712/mi is the 
mass ratio in the sense that mi is the more massive star, P is the 
orbital period and G is the gravitational constant. 

IT Com is less massive, the K giant primary is about 1.4 in- 
stead of 2.5 solar masses (TablefJJ, the secondary component could 
not be an F main sequence star like in 73 Leo (which would have 
about the same mass and would be visible in the optical spectra). 
From the rotational period and the radius of the primary (see Table 
[T]and TablefJ} the predicted rotational velocity is ~8 km s~ . The 
measured vsini is «5-6kms -1 , which implies a rotational incli- 
nation of about 40-45°. This, in turn, translates into a secondary 
mass of about 0.3 solar mass with the help of Fig. lATl i.e., the sec- 
ondary component is most likely an M dwarf. Taking into account 
the fairly high amplitude rotational modulation, the inclination an- 
gle can not be as low, so a coplanar orbit is probable, contrary to 
the result of Glebocki & Stawikowski dl997fc . 

B Masses of the secondary stars 

From the same calculations as above in appendix [A] the possible 
values of the secondary masses of V2253 Oph, IT Com and IS Vir 
can be found from Fig. IB II The estimated ranges of the primary 
masses and the orbital inclinations can be taken into account. The 
primary masses are as follows: V2253 Oph: 3 M®, (range 2.4— 
3.8 M ), IT Com: 1.4 M , (range 1.2-1.6 M ), IS Vir: 1.3 M , 
(range 1.1-1.5 M ). Eccentricities and radial velocity curve am- 
plitudes are given in Table [T] 
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Fig. Bl The relation between the orbital inclinations and 
mass ratios of V2253 Oph, IT Com and IS Vir. Solid lines 
mark the mean values of the primary masses, dotted lines 
indicate the mass ranges. 
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